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Abstract 

Two methods the complex energy shell model (CXSM) and the complex scaling (CS) approach 
were used for calculating isobaric analog resonances (IAR) in the Lane model. The IAR parameters 
calculated by the CXSM and the CS methods were checked against the parameters extracted from 
the direct numerical solution of the coupled channel Lane equations (CC). The agreement with the 
CC results was generally better than 1 keV for both methods and for each partial waves concerned. 
Similarities and differences of the CXSM and the CS methods are discussed. CXSM offers a direct 
way to study the configurations of the IAR wave function in contrast to the CS method. 

PACS numbers: 21.10.Sf,24.30.-v,24.10.Eq,21.60.Cs 
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I. INTRODUCTION 



The isobaric analogue resonances (IAR) have been discovered about four decades ago. 
They are the consequences of the approximate isospin symmetry. In the absence of isospin 
symmetry braking forces the analogue state would be degenerate with corresponding state 
of the parent nucleus. The Coulomb forces and isospin dependent interactions shift up the 
position of the analogue state and the analogue state may acquire a small width and may 
become a resonance. 

The IAR has attracted interest in the latest development of the study of the light exotic 
nuclei. The unusual properties of neutron reach nuclei can provide insights into the nuclear 
structure far from the valley of stability. The extreme neutron to proton ratios may help to 
understand the nuclear matter at extreme conditions. However the experimental study of the 
neutron rich nuclei around the neutron drip line are difficult. Since the IAR has essentially 
the same structure as the parent state it was suggested that instead of the neutron reach 



exotic nuclei ( n Li, 14 Be, 7 He, 9 He) 



their less exotic analogue states should be 



studied (with inverse kinematics) to gain information about the properties of these exotic 
nuclei. 

Recent developments in experimental facilities opened the possibility for identifying large 
number of exotic nuclei. To understand the structure of these nuclei new theoretical methods 
have been introduced for describing the dynamics of weakly bound or unbound nuclei from 
which nucleons can be emitted. Some of the new methods is e.g. the Shell Model in 
the complex energy plane (CXSM) 0, [f| or the Gamow Shell Model [7, 0, 10, ll| use 
the Berggren basis 12J. In the Berggren basis bound and resonant states are treated on 
equal footed and scattering states taken along a contour L of the complex energy sheet are 
included as it will be discussed later in detail. In the last few years this basis has been 



used successfully in a serious of works 
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17|. Since the extended use of this 



basis started not very long ago therefore we think that it is worthwhile to accumulate more 
experience concerning its accuracy and its dependence on their parameters. 

Another well established method for calculating resonances is the complex scaling (CS) 



method. CS has a strict mathematical foundation given in Refs. 



applications and the details of the CS method are reviewed in [21 
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The IAR states phenomenologically can be described by the Lane equations [23(] or they 
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can be studied by microscopically [2J|. Coupled channel (CC) Lane equations offer a simple 
but not trivial (the simplest multichannel) example on which both the CXSM and the CS 
approaches can be checked. 

Our aim in this work is to compare the parameters of the IAR calculated by different 
methods. We shall compare the characteristic features of the two methods working on the 
complex energy plane. The conclusion of such a methodical work can be useful later in 
analyzing experimental data in more realistic calculations. 

The CS can be applied only for dilatation analytic potentials and interactions. However 
some of the widely used potentials in nuclear physics are not dilatation analytic or dilatation 
analytic only in a limited range of the rotation angle, i.e. below the critical value of the 
angle. Therefore we repeated our calculation with a slightly modified Coulomb potential 
which is dilatation analytic. This comparison is very useful to compare the accuracy of the 
CXSM and CS methods if they give the same results in cases when both methods can be 
applied. 

In section III Al we summarize the features of the Lane equations. In section III Bl we 
describe the approximate solution of the Lane equations using the CXSM, while in section 
III CI we make a short description of the solution of the Lane equations using CS. In chapter 
III 1 1 we give the numerical results of the calculations. In the first part of chapter II III we 
compare the positions of the poles of the S-matrix calculated by the CXSM method with 
that extracted from the solution of the Lane equations. The results of the CS method are 
also presented here. The similarities and differences of the CXSM and CS methods are also 
discussed in that chapter. Finally in the last chapter we summarize the main conclusions of 
the paper. 



II. RESONANCE SOLUTION OF THE LANE EQUATION 

The Lane equations in the simplest case describe the quasi-elastic scattering of a proton 
and the IAR. We assume that the target nucleus has mass number A = N + Z and charge 
number Z. The ground state of the target has isospin T4 and isospin projection T 3 = = 
Ta- The target is bombarded by a beam of protons. 
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A. Lane equation 



The Hamiltonian of the target plus nucleoli system H can be divided to a part describing 
the internal motion of the target H(£) and their relative motion H re i 



H = #(£) + H rel . 



(1) 



The internal state of the ground state of the target is denoted by \A) and this state is the 
solution of the equation H(£)\A) = €a\A). The analog nucleus is denoted by A it has the 
same isospin Ta and isospin projection T^ — l. It is an excited state of the isobaric nucleus 
with Z + 1 protons and N — 1 neutrons. If we neglect the mass difference between the 
neutron and proton and denote the additional Coulomb energy of the analogue nucleus by 
A c then the eigenvalue of the internal motion of the analogue state is simply €a + A c and we 
have H(£)\A) = (e^ + A C )\A). Let \pA) and \nA) be the states formed by adding a proton 
and neutron to \A) and \A), respectively. The total wave function of the system may be 
written in the form 

* = \A)<f> p (r) + |i)0 n (r), (2) 



where <f> p (r) and 



(r) describe the relative motion. The relative motion part of the total 
Hamiltonian can be cast into the form 



rel 



K + V (r) + t • f Vi(r) + - t 3 )V c (r) 



(3) 



where K is the kinetic energy operator of the relative motion, Vq comes from the interactions 
independent form the isospin, Vc is the nuclear Coulomb potential and t • TVi(r) is the 
symmetry term due to the isospin dependent strong interactions. The vector operators t 
and T are the isospin operators of the nucleon and that of the target. Substituting the 
ansatz (j5J) into the Schrodinger equation H^> = S^f and taking into account the form (j3J) of 
the relative Hamiltonian we get the Lane equations [23 ]: 



K + Vo-^TjM + Vc-S, 



K + V Q + -{TA-VjVr - {£ p - A c ) 



+ \ ~T A V l( f> n = 



-T A V, 



(4) 



where £ p = £ — e A is the center-of-mass energy of the relative motion of the proton plus 
target system and the relative energy for the neutron plus analogue nucleus is £ p — A c . If we 



assume that the interactions are spherical symmetric then the relative motion can be divided 
into partial waves and there are no coupling between different partial waves characterized by 
orbital / and total angular momentum j quantum numbers. The numerical solution of the 
Lane equation has been carried out by using fourth order Runge-Kutta method. At each 
real £ p values we calculated two linearly independent solutions of the coupled equations. 
The physical solution with components cj) p and <p n being regular at r = was combined from 
these independent solutions. These components <j) p and <j) n were matched to the scattering 
(or outgoing wave) solutions of the corresponding channels at a distance where the nuclear 
potentials are cut to zero. This is a standard method described e.g. in Ref. [25]. 



B. CXSM a solution using Berggren basis 



In this section we calculate the complex energy eigenvalues of the IAR by diagonalizing 
the Hamiltonian ([T]) in combined Berggren bases 12j of the target plus proton and analog 
plus neutron systems. First we describe the Berggren basis for the protons. We consider 
an auxiliary problem, a radial Schrodinger equation with the diagonal potential of the first 
equation of the Lane equation (TJJ 



Kt + Vo- -T A V l + V c 



vg>{r) = 0, 



where 



Ki = - 



2fJL 



2 r ^2 



dr 2 



1(1 + 1) 



(5) 



(6) 



The discrete bound and resonance solution with energy En are denoted by Un (r) and 
the scattering solutions by E). Sometimes when it is obvious we will use the wave 

number k instead of the energy E and the scattering states along the contour in the lower 
half of the second energy sheet will be denoted by u^(r, k). 

The main advantage of the Berggren basis is that the single-particle basis set consists not 
only of bound states but also of poles of the single-particle Green function on the complex 
energy (wave number) planes and a continuum of scattering states taken along a complex 
contour L. A typical contour of the complex wave number plane is shown in Fig. [U The 
L + part of the contour goes from the origin to infinity in the lower half of the second energy 
sheet, while the L~ part of the contour makes exactly the same tour on the first energy 



sheet. It was observed in Refs. 
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271 ] that the contour L + need not return to the real 
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FIG. 1: Positions of the bound (61, 62) and decaying resonant (di, efo) poles of S(k) on the complex 
fc-plane and a possible choice of the complex contour L . 

axis at infinity. The shape of the chosen complex contour L = L + + L~ regulates which of 
the poles should be included into the Berggren basis forming the completeness relation of 
Berggren: 



5( r -r')= V u%\r)u%\r') + 2 / dku {p \r, k)u {p \r\ k) . (7) 

In this relation (and later) the notation n = b,d means that the sum over n runs through 
all bound states plus the decaying resonances lying between the real energy axis and the 
integration contour L + of Fig. HJ The integral in Eq. (JTj) is over the scattering states along 
L + (factor 2 is due to the symmetry of L). The poles denoted by d in the basis generally 
correspond to decaying resonances lying on the fourth quadrant of the complex fc-plane. 



The completeness relation in Eq. 
and it has been shown later in Ref. 



(171 ) was introduced for charge less particles in Ref. 



12 



28j that it is valid even for charged particles. Berggren 



completeness can be generalized by using a contour of different shape in which antibound 



states 29| lying on the negative part of the imaginary fc-axis are included in the sum in Eq. 
(JTJ) . Since the inclusion of antibound states is not optimal as far as the number of basis 
states is concerned [3o| we are not using antibound basis states in this work. 

Berggren introduced a generalized scalar product between functions defining a special 
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complex metric of the Berggren space 12[]. In the generalized scalar product in the left 
(bra) position of the scalar product the mirror partner state (denoted by tilde over the 
state) is used. This state corresponds to the reflection to the imaginary fc-axis. Due to 
this reflection in this scalar product in the integration the radial wave functions appears 
instead of the complex conjugate of the radial function. (This causes no difference for 
bound (and anti-bound) states lying on the imaginary axis.) This is the only modification 
in the scalar product since the spin-angular degrees of freedom remains unchanged. If the 



radial integral to be ca^ 



culated has no definite value then a regularization procedure has to 



be applied. Zel'dovich [3l| and also Romo [32|] suggested regularization methods but we use 



the complex rotation of the radial distance r beyond the range of nuclear forces [33|. 

The upper half of the complex /c-plane maps to the physical (or first) Riemann-sheet of 
the complex energy E ~ k 2 . The pole wave functions of this sheet are square integrable 
functions belonging to bound states. While the lower half of the complex fc-plane maps 
to the unphysical (or second) Riemann-sheet of the energy E. The pole wave functions of 
the second sheet are not square integrable functions and they belong to decaying/ capturing 
resonances lying on the lower/upper part of that energy sheet or antibound states lying 
on the negative real energy axis. The calculation of integrals in which these radial wave 
functions appear might need the use of regularization procedure. 

Since the number of basis states has to be finite the complex continuum has to be dis- 
cretized. It is preferable to use as discretization points E^ the abscissas of a Gaussian 
quadrature procedure and the corresponding weights of that procedure are denoted by hi. 
By discretizing the integral in Eq. (JTj) one obtains an approximate completeness relation for 
the finite number of basis states: 

M 

S ( r - r >) « ^1%, E^)w^(r>, (8) 

n=b,d,c 

where c labels the discretized contour L + states. If E$ corresponds to scattering en- 
ergy from the contour L + then the scattering state of the discretized continuum is denoted 
by Wn\r, E&) = VKv£(r, En^) and if E^ corresponds to a normalized pole state then 
Wn\r,En^) = Un\r). The set of Berggren vectors form a bi-orthonormal basis in the 
truncated space 

<tB«|«;W >=*„,„. (9) 



7 



The Berggren basis for neutrons is defined similarly but the auxiliary problem uses the 
diagonal part of the second equation of (j3J) 



K l + V Q + ^{T A - \)V x -EW 



«W(r) 



0. 



(10) 



Having fixed the Berggren basis for neutrons and protons we take the ansatz ([2]) but the 
relative motion functions are expanded on the corresponding Berggren basis 



M v 



and 



_%=\ 

' Mn 

i=i 



y 



Ijm 



(n)^ 



y 



Ijm > 



(11) 



(12) 



where 3^zjm denotes the spin-angular part of the wave function. Using Eq. (j3j) we get the 
following set of linear equations for the unknown complex expansion coefficients C 4 and 



a 



(») 



« J - W + ££i<tf?V«WW = o 



,W\^(n) 



(13) 



and 



(e^ - (s p - A c ))cr + E7=i^r\H^')c^ = o 



t(n) 



,T.(")| 



,(p)\rr(p) 



A- 



1, . . . , M n , 



(14) 



where the coupling potential is 5v = y^Vy. The above two equations can be combined 
into one matrix eigenvalue equation with dimension M p + M n . By diagonalizing the matrix 
of the Hamiltonian we get M p + M n complex eigenvalues 8 p v = 1, . . . M p + M n . One of the 
complex eigenvalues is identified by the energy of the IAR. The identification in general 
is easy because most of the other unbound states correspond to discretized contour states 
and they are lying far from the position of the IAR at Sjar 
function of the IAR the dominant component is a bound neutron state. 



E r — ztJ and in the wave 



C. Complex scaled Lane equation 



The poles of the Green-operator on the complex energy plane can be determined with the 
help of the complex scaling. The CS mathematically well founded 18|, [l9|, l20| and has many 
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applications in atomic, molecular and nuclear physics. We demonstrate the effect of the CS 
on an example of a single-particle Hamiltonian h. The real angle 9 of the CS rotates the 
coordinates of the particle to complex, i.e. r is simply replaced by exp(i9)r. More precisely 
the effect of the CS can be given with the help of an operator U(9). It acts on an arbitrary 
function g(r) as 

U(9)g(v) = ex P (i~%(re* ) . (15) 
The complex scaled Hamiltonian is the following 

h e = UehU,- 1 . (16) 

The kinetic energy K = — f^A r transforms due to the complex scaling to 

U e kU d - 1 = exp(-i20)(-|U r ) , (17) 

and a local potential V"(r) transforms to the form: 

V 9 (r) = UeVir)^ 1 = V(rexp(i9)) . (18) 

Assume that x^( r ) is a bound or resonance eigenfunction of the the Hamiltonian h and 
the corresponding eigenvalue is E u then the function xti r ) = exp(i|^)x I/ (re* el ) will be the 
eigenfunction of the complex scaled Hamiltonian h$ with the same eigenvalue E v . The 
advantage of the CS is that the function xii r ) is square integrable even if the original state 
was a resonance wave function. The square integrability of xii r ) allows that it can be 
approximated well with finite expansion using only square integrable basis functions. 

The Lane equation can be considered as an eigenvalue problem of a two by two matrix 
Hamiltonian 

\ K + V + ^V 1 + A C 

The Lane-equation (j3J) can be cast into the form 

/ A, (*\ \ 

(20) 





The generalization of the operator Ug is straightforward 



and the complex scaled matrix Hamiltonian is Tie = Ug7iU e 1 . The eigenvalue problem of 
this operator 

(22) 











He \ 




Ml 






V <{*) ) 


V ) 



in components gives the following set of equations 



[H e p -^ P ]< + \^ A VU e n = 



(23) 



where H° = exp(-i29)K + V 9 - \T A Vf + V% and H e n = exp(-i29)K + V 6 * + \{T A -l)Vf ' + A c . 
We will refer to (123]) as complex scaled Lane-equation. Since the functions 4> p {y) and <j> n (r) 
are square integrable we can make the approximation 



E^f* 



and 



Ijm 



i=l 



3>1 



Ijmi 



(24) 



(25) 



where V^tV) and V'j ( r ) are arbitrary square integrable basis functions. Substituting these 
forms into the (123|) we get a matrix eigenvalue equation. In details we have 



M n 



m=l 



m=l 

= ^E<A W I^ ) >C^ A; = 1,...,M P , 



(26) 



m=l 



and 



m=l 



m=l 



(27) 



m=l 



The solution of these equations provide us M p + M n number of complex eigenvalues. The 
majority of these eigenvalues correspond to the discretization of the rotated down continuous 
spectrum. The bound and resonance poles can be clearly identified and the accurate value 
can be determined using the so called 6 trajectory technique. 
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III. NUMERICAL RESULTS 



We applied the methods described in section iHBl and HTCl for the description of the IAR- 
es in the 209 Bi nucleus with large neutron excess. We studied several analogue resonances 
in the p + 208 Pb system. For illustrative purposes we selected an IAR which in our simple 
model is the analog of the ground state of 209 Pb, i.e. a g$/2 single particle state. The effect of 
the double magic core is described by a phenomenological potential. We used Woods-Saxon 
(WS) forms for both the diagonal and the coupling potentials in (j4]). The WS forms cut to 
zero at a finite distance: R ma x — 20 fm 

f V ws (r) if r < Rmax 
V™ s (r) = 1 K ' max (28) 

^ if r > R max . 

The spin-orbit part of the potential has the usual derivative form: 

V™(r) = 2(f- s) 6 ^ R . (29) 

ra 7 (i + e ^)2 

It is also cut to zerojit R ma x- The numerical values of the potential parameters were taken 
from an early work 



34j. For the sake of simplicity the radii and the diffuseness were taken 
the same values for protons and neutrons and for the common spin-orbit term: r = 1.19 
fm and a = 0.75 fm, V s0 = 11.6 MeV. For Coulomb potential we assumed that the charge 
Ze of the target is homogeneously distributed inside a sphere with radius R c = rcA 1 ' 3 with 
sharp edge 

, f oS-[3- (ir) 2 ] if r < R c 
V c (r) = Ze 2 { 2Rj KRc)l (30) 

I 1/r if r > R c . 

The depth of the nucleon potential was 56.4 MeV and the strength of the symmetry potential 
was 0.5 MeV. Therefore the diagonal WS potential felt by the proton was 61.9 MeV and by 
the neutron 51.15 MeV according to the Lane equations in Eq. The Coulomb radius 
was identical with the one of the nuclear potential. The Coulomb energy difference was also 
the same as in Ref. Q A c = 18.9 MeV. 

In the CXSM the elements of the single particle bases are calculated in the diagonal 
potentials appearing in the corresponding channels of the Lane equations. The single particle 
energies for the gg/2 neutron and proton orbits are summarized in Table [B The vertexes 
of three different proton and neutron L contours for the gg/2 case are shown in Table HT1 
The numbers of the discretization points N{ of the segment [V£, Vi+i] are shown between the 
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State 


neutron 


proton(Eq. <[5D])) 


proton(Eq. (1321)') 


lff9/2 
2 59/2 


(-22.878,0.0000) 
(-4.060,0.0000) 


(-11.894,0.0000) 
(7.674, -6.210- 4 ) 


(-13.975,0.0000) 
(6.070, -2.010- 5 ) 



TABLE I: Energies of the discrete g 9/ / 2 single-particle basis states for neutrons and for protons 
corresponding to the Woods-Saxon potential described in the text. The Coulomb potential for 
protons is either the usual one in Eq. (|30p or the dilatation analytic one in Eq. (|32p . Energies are 
in MeV. 

vertex points. To calculate IAR-es we could use neutron contours taken along the real axis. 
For other partial waves we used large variety of contours. 

In the CC method we solved the coupled Lane equations for a fine equidistant mesh of 
the bombarding proton energy £ p = S ,S + d£, ...,£ max in the center of mass system and 
calculated the scattering matrix elements for each energy values: S(S P ). 

In order to determine the parameters of the IAR we fitted the tabulated values of the 
S(Si) by the following form: 

S(S P ) = e 2i5 ^\l - i—^ ) , (31) 

£p — £>IAK 

where Sjatz denotes the complex energy of the IAR, i.e. £iar = E r — i^, T is the full width 
and T p is the proton partial width of the IAR. Below the threshold of the 208 Pb(p,n) 208 Bi 
reaction the total width is equal to the the partial width: r = T p in our model. Equation 
fl3Tj) represents a one pole approximation to the S-matrix in the proton channel. 

For the background phase shift in the entrance channel S p (S p ) we take a linear energy 
dependence in order to better reproduce the non-resonant background. Naturally all these 
quantities refer to definite l,j partial waves. The best fit parameter values are listed in the 
E r (CC) and T(CC) columns in Table HLT1 The one pole formula of (}3~Tj) gave excellent fit to 
the tabulated values of S(£i) in all cases in Table IIHI 

The numerical values of £iar are shown as E r (C X S M) and T(CXSM) in Table IIHI 
in comparison with the results extracted from the solution of the coupled Lane equations 
E r (CC) and T(CC). As one can see from the comparison the positions and the widths 
calculated by the CXSM agree well (within 1 keV) with the result of the CC. 

Let us discuss below briefly how this agreement has been achieved. We optimized the 
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proton 


neutron 
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n 
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n 
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TO 
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nn 
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v 2 


^OU,-U.4L J 


^ou,-z.uyo j 


("in n\ 


{OU,U) 




Cm ~\(\\ 




o 
Z 


4 


o 
z 


n 



n 






v 3 


(30,0) 


(100,-6.993) 


(100,0) 


(100,0) 




(10,0) 


N 3 







4 








4 


v 4 


(100,0) 




(200,0) 


(200,0) 




(30,0) 




















v 5 


(200,0) 










(100,0) 



TABLE II: Integration contours for g 9 / 2 protons and neutrons given by vertexes Vi (in MeV) and 
the number of Gaussian points iVj. The iVj values are the ones necessary to reach the 1 keV 
accuracy for the IAR or for the broad resonance at £ p = (23.996, —6.147) MeV. Contours LP3 and 
LN3 were used for the broad resonance. 

shape of the contours and the number of points along the contours separately for neutrons 
and protons and for different partial waves. The shape of the contour is fixed by the vertexes 
which were chosen to be able to include the narrowest single-particle resonant states. We 
observed that the contour should not go close neither to the energies of the resonances 
included into the basis nor the IAR resulted by the diagonalization. The last vertex point 
i.e. the energy of the last segment with iVj ^ was crucial to get good agreement for both 
the real and the imaginary parts of the IAR energy calculated by solving the Lane-equation. 

We tested the convergence of the IAR energies by increasing the number of discretization 
points and stopped to increase it when the energy did not change. After that we continued 
with the next interval and increased the points of that interval similarly. After going through 
all the intervals we optimized the number of mesh-points by reducing them until the energy 
in keV did not change. We also tested the convergence of the IAR energies by varying the 
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1 3 


E r {CXSM) 


E r (CC) 


T(CXSM) 


T{CC) 


99/2 


14.954 


14.954 


0.046 


0.047 


^11/2 


15.526 


15.526 


0.003 


0.003 


d-5/2 


16.445 


16.444 


0.141 


0.140 


Sl/2 


16.918 


16.917 


0.156 


0.156 


97/2 


17.367 


17.367 


0.086 


0.084 


ds/2 


17.441 


17.440 


0.144 


0.145 


J15/2 


18.774 


18.774 


0.006 


0.006 



TABLE III: Comparison of the IAR parameters calculated by using the CC and CXSM methods 
for different partial waves with Coulomb potential Eq. (|30p . Energies are in MeV units. 

positions of the vertexes. If the contour goes very far from the real axis i.e. if we choose the 
value of the imaginary parts of the vertexes considerably larger than the ones in Table HT1 then 
the degree of agreement might be spoiled even if we choose larger number of discretization 
points. We found for all partial waves that the IAR resonances are not very much sensitive 
to the low energy part of the continuum (below 5 MeV), neither for neutron nor for proton. 
At high energy however a cutoff smaller than 30 MeV affects the convergence of the IAR 
energy. 

In order to be able to compare pole solutions for calculation of the resonance parameters 
of the IAR we repeated the calculation by applying the CS for the solution of the Lane 
equations. Unfortunately the Coulomb potential of a charged sphere with sharp edge is 
not dilatation analytic because this form becomes discontinuous for 6^0. We used the 
Coulomb potential expressed by the error function which is dilatation analytic. This form 
of the Coulomb potential 

Vo(r) = Ze> Erf(r/a) (32) 



is widely used in both atomic and nuclear physics [26|, [35|. In the resonating group model it 
can be obtained as the direct folding interaction between nuclei [36|. The numerical value 
of the parameter a = 0.31 fm was adjusted to the Coulomb potential in Eq. fl30|) . For the 
nuclear potential we kept the WS form which is dilatation analytic until the rotation angle 
is below the critical angle: 6 < 9 crit = arctg(^). 
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For the solution of the complex scaled Lane equation we used the Laguerre mesh basis 
functions 

4 U \r) = i-m 1/2r -y^ exp(-r/2), (33) 

where v — p,n. The mesh points given by LuX r i) = 0; where Lm v {t) is the Laguerre 
polynomial. The advantage of this basis is that the matrix elements of any local potential 
is extremely simple (37)]. This type of basis functions are proved to be very accurate both 
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38 



391 ] . One can introduce an 



in simple model calculations and in three body problems 

additional simple scaling parameter of the basis j^] for this parameter we used the 0.3 fm 
value. 

The agreement between the pole positions calculated by the CXSM and the CS method 
is extremely good for all partial waves in Table IIV1 One can see in this Table that the 
agreement with the numerically exact solution of the Lane-equation (CC) is as good as in 
the previous case when the CXSM method was used with the standard Coulomb potential. 
The maximal difference does not exceeds 1 keV. 



unperturbed n states 
unperturbed p states 
perurbed states 



"'-20 -10 10 20 30 40 50 60 

Re(E) [MeV] 

FIG. 2: Positions of the g§/2 Berggren basis states (circles for neutrons and squares for protons) and 
the results of the CXSM method (filled circles) on the complex i?-plane for the Coulomb potential 
Eq. p}. 

To understand better the formation of the IAR let us consider again the I = 4, j = 9/2 
case as an example. In Fig. [2] we show the positions of the unperturbed states forming 
the Berggren basis (denoted by circles for neutrons and squares for protons) and some of 
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1 3 


E r {CXSM) 


E r (CS) 


E r (CC) 


T(CXSM) 


T{CS) 


T(CC) 


99/2 


14.933 


14.933 


14.934 


0.041 


0.041 


0.042 


^11/2 


15.493 


15.493 


15.493 


0.002 


0.002 


0.002 


d-5/2 


16.436 


16.436 


16.444 


0.121 


0.120 


0.120 


Sl/2 


16.913 


16.913 


16.913 


0.127 


0.127 


0.128 


97/2 


17.350 


17.349 


17.349 


0.076 


0.075 


0.074 


ds/2 


17.434 


17.434 


17.433 


0.118 


0.119 


0.120 


J15/2 


18.752 


18.751 


18.752 


0.005 


0.005 


0.005 



TABLE IV: Comparison of the IAR parameters calculated by using the CC, CXSM and CS methods 
for different partial waves with the dilatation analytic Coulomb potential Eq. (|32p . Energies are 
in MeV units. 

the results of the diagonalization (perturbed states denoted by filled circles) on the complex 
energy plane. A section of the real £"-axis and of the lower half of the complex plane is 
shown. In this case the neutron contour is along the real axis while the proton contour has 
a trapezoidal shape with vertexes denoted by LP1 in Table [Til In order to see better the 
region of our interest the states with energies higher than 60 MeV are not shown in FigJ21 
The discrete basis states are listed in Table [TJ The bound basis states are the lgg/2 proton 
state at E^' = —13.975 MeV and the two neutron states: lgg/2, 2g 9/ / 2 at E\ n ' = —22.878 
MeV and at E\ n) = -4.060 MeV which are shifted up by A c = 18.9 MeV . The narrow 2g 9/2 
proton resonance at E± = (6.070, —2 x 10 -5 ) MeV seems to lie on the real axis. 

Most of the perturbed states lie close to the positions as the corresponding basis states 
since the coupling symmetry potential term cases only a small shift for these states. One 
of the exception is the IAR at Sjar = (14.933, —0.021) MeV which shifted down well below 
the bound 2gg/ 2 neutron state which is the main component of its wave function C^ 1 ' = 
(0.9921, —0.0047). The second largest component is that of the 2g 9 / 2 proton resonance with 
= (—0.1194,0.0002). The other perturbed states which do not fit to the path of the 
contours are states based on contour states but fall off the contour because of the finite 
number of discretization points. If the number of discretization points are increased they 
move closer to the contour. They move also with the contour if we change the shape of 
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the contour in contrast to the IAR which remains in the same position. Of course the IAR 
should lie above the proton contour in order to be explored. This feature is very similar to 
the one observed in the CS calculation. 

From the mathematical theory of the complex scaling 
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22] it is known 



that the continuous part of the spectrum of the complex scaled Hamilton operator consist 
of half lines on the complex energy plane. The half lines start at the thresholds and they 
are rotated down from the real axis by 29. In our calculation we have used M p = M n = 
100 basis functions and received two hundred approximate complex eigenvalues from the 
diagonalization. These eigenvalues are plotted on Fig. [3] for two different 9 values 9 = 2° and 
9 = 4°. From this Figure it is obvious that the vast majority of the eigenvalues correspond 
to discretization of the continuous spectrum. However there are a few eigenvalues which are 
independent from the complex scaling parameter 9. These are denoted by letters 61,62 and 
ri,r 2 in the Figure [3j The r 2 is the IAR which is based mainly on the 2 g 9 / 2 bound neutron 
state as we have seen in the CXSM calculation before. The r± resonance is based mainly 
on the narrow 2g$/2 proton resonance at = (6.070, -2 x 10~ 5 ) MeV. The bound states 
bi and 62 originate on the I5V2 proton state at J5 f = —13.975 MeV and the lgg/2 neutron 
state at E\ n) = -22.878 MeV which are shifted up by A c = 18.9 MeV. 



The lower part of the Figure [3] show the so called 9 trajectory i.e. the complex energy 
plane in the vicinity of the IAR when the complex scaling parameter changes between 9 = 2 
and 9 = 8 degrees with step size of one degree. There is a small change in the position 
and width of the resonance (this should be independent form the value of 9) but this comes 
from the fact that a finite basis is used. This phenomena is well know in all complex scaling 
calculation and there are methods how to select the best approximation for the resonance 
{40! . The resonance position and width values given in Table IIVI correspond to calculations 
with 9 = 4°. 

A similarity of the CXSM and the CS method is that results become less accurate if the 
contour of the CXSM or the rotated half lines are lying close to the resonance. To get high 
accuracy the resonance has to be well explored i.e. should lie far above the contour. The 
rotated half lines of the CS play similar role as the contours of the CXSM therefore we shall 
call the half lines of the CS method also contours. Only the resonances above the contours 
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FIG. 3: Upper part: positions of the states on complex .E-plane with CS method with the 
rotation angle 9 = 2° (filed circles) and with 9 = 4° ( open circles) for the Coulomb potential Eq. 
(|32p . The IAR is the state r-i- Lower part: the vicinity of the IAR, the CS scaling parameter 9 is 
varied between 2° and 8° with a step size of 1°. 

can be calculated. This means that the 3gg/2 neutron resonance at = (4.929, -6.035) 
ar g(£'( n )) = 50.76° or the corresponding perturbed solution can not be calculated by the 
CS method since they can not be explored due to the critical angle of the WS potential 
Qcrit — 18.48°. It can be calculated however by the CXSM using contours LP3 — LN3 and 
we get for the perturbed energy £ p = (23.996, —6.147) MeV arg(£ p — A c ) = 50.33°. 

The similarity of the methods can be seen even better if we try to use a contour in 
CXSM, which resembles to the rotated continuum of the CS calculation. In Fig. H] we 
present the results of the CXSM calculation in which the contours LP2 and LN2 were 
chosen to be the same as the one corresponding to the optimal 9 = 4° rotational angle of the 
CS calculation. One can see that the IAR is well separated from the two contours starting 
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at the origin and the one starting at the neutron emission threshold. The unperturbed pole 
closest to the IAR is the bound 2g 9 / 2 neutron state which is the dominating component 
of the IAR wave function with amplitude: C\ n ^ = (0.9918,-0.0043). The second largest 
component of the IAR wave function is the one of the 2gg/ 2 proton resonance with amplitude 
C\ p) = (-0.1194,0.0001). The wave function of the IAR practically unchanged as far as the 
discrete components are concerned with respect to the case with the contours used in Fig. 
M (LP1 and a real neutron contour). The energy of the IAR is S p = (14.93309, -0.02058) 
MeV coincide with the one S p = (14.93303, —0.02062) MeV with the contours used in Fig. 
|2] within the numerical error 1 keV estimated from the deviation from the CC results in 
Table IIVI This good agreement convinces us that the use of the LP2 and LN2 contours 
which resembles to the contour of the CS could also be used for calculating the IAR. The 
components of the different scattering states taken from the different contours are certainly 
very different but the summed contribution of the proton and neutron contours are basically 
the same. Since both are small numbers their numerical values have little importance. For 
the gg/2 IAR the neutron continuum has negligible effect. For other partial waves this effect 
is also small but not completely negligible. 

An important difference between the results presented in Figs. [3] and H] is that in Fig. 
[3] only perturbed states are shown since in the CS method the basis states used are not 
eigenstates of any unperturbed Hamiltonians. Therefore from the coefficients of the wave 
function of the IAR explored we can not estimate the role of the unperturbed neutron and 
proton states easily. To get similar quantities we have to calculate the unperturbed state 
with the same CS contour and we have to calculate overlaps with the IAR wave function. 

IV. SUMMARY 

Let us summarize briefly the results we received in this study. We reproduced the results 
of the direct numerical solution of the coupled Lane equations by diagonalizing the Hamil- 
tonian in the full n-p Berggren basis i.e. using the CXSM method. The IAR parameters 
were extracted from the S(£ p ) calculated by solving the Lane equation along the real £ p -axis 
by fitting it using the one pole approximation Eq. ( J3T1) . The fitted position E r (CC) and 
the width T(CC) of the IAR was compared to the result of the CXSM calculation and the 
agreement was generally better than 1 keV for all partial waves in which we had IAR. In 
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FIG. 4: Positions of the Berggren basis states (circles for neutrons and squares for protons) and 
the results of the CXSM method (filled circles) on the complex i?-plane for the Coulomb potential 
Eq. (|32p with contours LP2 and LN2 in Table [III They are similar to the optimal contour of the 
CS method. 

the wave function of the IAR furnished by the CXSM the contribution of the bound neu- 
tron state has the dominant role and the proton resonance has a non negligible effect. The 
integrated effect of the proton continuum is small but essential to produce the correct width 
for the resonance. We studied the details of the different parts of the continuum segments 
and the necessary numbers of the discretization points on the different segments. The role 
of the cut off energy and the low energy part of the continuum were also investigated. The 
neutron continuum played very small effect for the IAR-es. 

The pole position of the IAR was calculated by complex scaling method as well. For 
that we modified the Coulomb potential for a dilatational analytic one and repeated the CC 
calculation and CXSM method with the modified Coulomb potential. We received very good 
agreement to the numerical solution of the coupled Lane-equations both with the CXSM 
and the CS methods. Therefore we conclude that in this case the CXSM and the CS method 
give basically the same results apart some numerical errors which naturally not the same 
in the two type of calculations. This agreement suggests that the two method are basically 
equivalent in those cases when both methods can be applied. 

Besides the similarities and differences of the CXSM and the CS methods discussed so 
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far there are further important differences between them. The application of the uniform 
CS method used here is restricted to dilatation analytic potentials and the range of the 
rotational angle could also be limited. On the other hand in the CXSM method the shape 
of the contour can be chosen with large flexibility although to go too deep into the complex 
energy might spoil a bit the accuracy of the calculated results. Another advantage of the 
CXSM is that the structure of the resonant state can be seen directly from the coefficients 
of the perturbed wave function. In the CS method the same information can be explored in 
a more indirect way. 
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